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ABSTRACT 

We discuss properties of oscillatory convective modes in low-mass red giants, and 
compare them with observed properties of the long secondary periods (LSPs) of semi¬ 
regular red giant variables. Oscillatory convective modes are very nonadiabatic g~ 
modes and they are present in luminous stars, such as red giants with logL/pQ > 3. 
Finite amplitudes for these modes are confined to the outermost nonadiabatic layers, 
where the radiative energy flux is more important than the convective energy flux. 
The periods of oscillatory convection modes increase with luminosity, and the growth 
times are comparable to the oscillation periods. The LSPs of red giants in the Large 
Magellanic Cloud (LMC) are observed to lie on a distinct period-luminosity sequence 
called sequence D. This sequence D period-luminosity relation is roughly consistent 
with the predictions for dipole oscillatory convective modes in AGB models if we 
adopt a mixing length of 1.2 pressure scale height {a = 1.2). However, the effective 
temperature of the red-giant sequence of the LMC is consistent to models with a = 1.9, 
which predict periods too short by a factor of two. 

Key words: convection - stars: AGB and post-AGB - stars: late-type - stars: low 
mass - stars: oscillations . 


1 INTRODUCTION 


Long period variables (LPVs) are red giant stars which 
are known to obey several period-luminosi t y relation s (e.g. , 


Wood_eyaJj_ 199S^ Kiss fc Beddind l2003l : I It a et all l2004l : 
Soszvriski et al. 20041 . 20071) . Among these relations, the so- 


minosi ty 

[iooi; 11 


called sequence D relation corresponds to the longest pe¬ 
riod for a given luminosity, with periods ranging from about 
200 d to 1,500 d. This period is actually a long secondary pe¬ 
riod (LSP) co-existing with a primary period which is about 
8 times shorter an d which typically falls on sequence B (e.g., 
IWood et al.lll99^ ). Although other sequences (except for se¬ 
quence E) a re generally associa t ed with radial or nonradial 
pulsations (IWood et al.l 19991: iDziembowski fc Soszvhskil 


201C ; iTakavama. Saio &: Ital 20131: iMosser et al.ll2013l: IWoodI 


20151 ). the origin of sequence D is unknown . Sinc e 
the discovery of sequence D by IWood et all lll999l) . 
its properties and origin have been discussed widely 
(e.g.. Wood. Olivier fc Kawalerj|2004l : ISoszvhski et ahlkoOTl : 


iNicholls et all l2009l ). Various models have been consid¬ 
ered such as semi-d etached binaries dWood et al.l 1 19991 : 
[Sosz^nski_e£_^ 2007 ), rotatin g spots with dust formation 
llTakava ma^_Wood_&^a 20151) . and radial and nonradial 
pulsation dWood et all 1 19991) . However, none of them give 
a consistent explanation for the LSP variations. 


There is considerable evidence that radial pulsation 
is not the main mechanism associated with the LSPs. 
The observed radial velocity variations of a few k m s“^ 
(|Wood. Olivier fc Kawal^l2004l : iNichoIls et al.l[2009l ). com¬ 
bined with the long periods, would produce a large frac¬ 
tional radius change if the LSP is caused by radial pulsa¬ 
tion. However, there is no direct evidence, for example from 
changes in effective t emperature, for a large radius change 
dNicholls et al.l [20091) . In addition, no change in the pri¬ 
mary period of the star is observed as might be expected 
if th e radius is changing significantl y with the LSP varia¬ 
tion dWood. Olivier fc Kawal^l2004l ). Finally, we note that 
the fundamental mode radial p ulsation period of models is 
much shorter than the LSP (e.g. I Wood et al.l 19991 ). All these 
factors indicate that the LSPs are not caused by radial pul¬ 
sation. 


One previous s uggestion has a ssociated convection with 
LSPs in red giants. IStother j (l2010l) considered giant convec¬ 
tion cells for the origin of the LSPs. He argued that the 
turnover time of giant convection cells is comparable with 
observed length of LSPs, and that the observed radial veloc¬ 
ity changes can be interpreted as convection motion, which 
causes no radius change. However, it is not clear why the 
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turnover time would be related to a distinct periodicity with 
a light amplitude of up to one magnitude. 

In this paper we discuss oscillatory convective modes in 
red giants. They are oscillations confined to the outermost 
layers of the deep envelope convection. Dipole oscillatory 
convective modes have properties favorable for the cause of 
the LSPs, such as periods longer than the fundamental radial 
mode and small temperature variations. We also note that 
the above mentioned difficulties in radial pulsations associ¬ 
ated with the observed radial velocity variations of LSPs are 
not applicable for a dipole mode pulsation because it does 
not change mean stellar radius. 

We discuss evolutionary models on which we performed 
pulsation analysis in Section [2] The properties of the oscilla¬ 
tory convective modes are discussed in Section [S] We make 
comparisons with observations of LSP variability in Section 

H 


2 MODELS 

Fig.m is a HR diagram showing the positions of sequence 
D stars (small dots), and some evolutionary tracks on the 
red-giant branch and asymptotic giant branch (AGB) for 
initial masses of 1.0, 1.3, and 2.0 M© with a = 1.9 and 
1.2, where a, a mixing-length parameter, is defined as 
the mixing-length divided by the pressure scale-height. Se- 
quence D stars are ado pted from OGLE III identification by 
dSoszvnski et al.ll2009ll . The luminosity and effective tem¬ 
peratures were obtained from 2 MASS photometry data of 
J and K llSkrutskie et al.ll2006ll and the temperature and 
bolom etric correction - color relations by Houdashelt et ahl 
(l2000ll and [Houdashelt. Bell fc SweigartI 1 200C ), an d a dis¬ 
tance modulus of 18.54 mag dKeller fc WoodI 2006l l of the 
Large Magellanic Cloud (LMC). The 2MASS J,K magni¬ 
tudes were converted to the GIT system using the relation 
derived bv ICarpenteil ll200ll L 

Stellar evolution models were obtained using the MESA 
(Modules fo r Experiments in S tellar Astrophysics; version 
7184) code JPaxton et al.l[2013h . The hydrogen and metal 
abundances adopted are {X,Z) = (0.73,0.008). Models of 
various initial masses ranging from 0.9 M© to 2.0 M© were 
computed with several values (between 1.0 and 1.9) of the 
mixing-length parameter a. 

The central helium hashes start to occur near the tip 
of the red giant branch (RGB) at logL/LQ « 3.3. Af¬ 
ter a few shell hashes, the star transits to a steady cen¬ 
tral helium burning around logL/L© 2. After cen¬ 
tral helium exhaustion, the evolution along the asymp¬ 
totic giant branch (AGB) starts. Here helium shell hashes 
(thermal pulses) occur in which the luminosity changes 
cyclically while the average luminosity increases gradu¬ 
ally. The evolution calculations were stopped when each 
model reaches at logL/L(7) ~ 4.0. Wind mass lo ss rate of 
de Jager. Nieuwenhuiizen fc van der HuchtI (Il988l l was in¬ 
cluded in the calculations. 

The majority of the sequence D stars are AGB stars, 
whose loci on the HR diagram look consistent with evolution 
models with a = 1.9. 





Figure 1. Selected evolutionary tracks computed by the MESA 
code for initial masses of 1.0 Mq, 1.3 Mq and 2.0 Mq with the 
mixing-length parameter, a = 1.9 and 1.2. The evolutions were 
started at a pre-main sequence model and calculated through the 
thermal pulsing AGB stage up t o logL/Lf?^ ^ 4.0. Sm all dots 
are sequence D stars identified bv ISoszvnski et al.1 ll2009ll (OGLE 
III). 


3 OSCILLATORY CONVECTIVE MODES IN 

RED GIANTS 

IShibahashi fc Osakil (Il98ll f discovered that nonradial g“ 
modes, whose frequencies are purely imaginary (meaning 
monotonic growth or decay) in the adiabatic condition, 
corresponding to the convective instability, become oscil¬ 
latory in extremely nonadiabatic condition. We call them 
oscillatory convective modes. They are highly unstable (or 
stron gly excited) with growth tim e comparable with the pe¬ 
riod. IShibahashi fc Osakil lll98lf) discovered the oscillatory 
convective modes in studying nonradial pulsations of high 
angular degrees ^ > 10 in luminous (10® Lq) models hotter 
than the cepheid instability strip. Such high degree modes 
are not expected to be observa ble because of cancellation on 
the surface. Thirty years later. [Sai3 (1201 ill found low degree 
£ < 2 examples of these oscillatory convective modes in hot 
massive stars where convection zones are associated with the 
Fe opacity peak at T ~ 2 x 10® K. 

In this paper we consider oscillatory convective modes 
in luminous red giants having deep convection envelopes. 
Frequencies and amplitude distributions in the stellar inte¬ 
rior w ere obtained by the method described in ISaio fc CoxI 
lll980li . in which the Lagrangian perturbation for the diver¬ 
gence of convective flux is neglected (i.e., SX-F^onv = 0). The 
temporal variation is expressed as exp(icrt) with complex an¬ 
gular frequency cr, so that an oscillation mode is excited (or 
overstable) if Ui <0, where Ui is the imaginary part of cr. 
To represent the value of eigenfrequency, we use normalized 
(non-dimensional) eigenfrequency cu defined as 
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Figure 2. Top panel: Amplitude distributions of typical oscil¬ 
latory convective modes of -^ = 1 (black lines) and £ = 2 (red lines) 
are shown as a function of log T in a model on the AGB branch 
for a initial mass of 2.0 Mq with a mixing-length parameter of 
a = 1.2. Solid and dashed lines are, respectively, the real and 
imaginary parts of the normalized radial displacement Sr/r. The 
real and imaginary parts of the normalized angular frequency lj 
and the period IT in days for each mode are shown. We adopt in 
this paper the convention that a pulsation mode is excited if the 
imaginary part is negative (i.e., < 0). Bottom panel: Some 

physical variables in the AGB model are shown, where rconv and 
Tth are, respectively, convective turn-over time and thermal time 
in units of days. 


UJ = 



( 1 ) 


where R, M, and G are stellar radius and mass, and grav¬ 
itational constant, respectively i.e., the pulsation frequency 
is normalized by the Kepler frequency at the stellar surface. 

The top panel of Fig. [2] shows the amplitude of the ra¬ 
dial displacement 5r/r for oscillatory convective modes of 
1=1 (black lines) and 2 (red lines) in the interior of an 
AGB model of an initial mass of 2.0 M© at logL/LQ = 3.5 
with a = 1.2. Solid and dashed lines are, respectively, real 
and imaginary parts of the radial displacement. Actual ra¬ 
dial displacement is given by the real part of [6, (f>) 

where YJ^{6, cp) is a spherical harmonic and 9 and (p are the 
co-latitude and azimuthal angle, respectively. 

The amplitude distributions are remarkably similar for 
the i = 1 and 2 modes. Finite amplitudes are mostly con¬ 
fined to the outer layers where the superadiabatic temper¬ 
ature gradient V — Vad is large and the thermal time Tth 
is shorter than the periods, indicating that the oscillatory 
convective modes are very nonadiabatic. Here, V and Vad 
are defined as 


„ _ d\nT ^ _ fd\nT\ 

^ dlnP’ ^ VainpJs’ 


( 2 ) 


Figure 3. Real (upper panel) and imaginary (lower panel) of 
normalized eigenvalue of oscillatory convective modes in AGB 
(non-flashing) models of each mass (black lines — 1.0 Mq; red 
lines — 1.3 Mq; blue lines - 2.0 Mq). Solid and dashed lines are 
for models with a = 1.2 and 1.9, respectively. 


where P and S are the pressure and the specific entropy, 
respectively. The thermal time rth is defined as 

(M - Mr)C,T 

Rh — j- , (3) 

where Mr is the mass within the sphere of radius r (i.e., 
M — Mr is the mass above the sphere of radius r), Cp is 
the specific heat at constant pressure, and Lr is the local 
luminosity at radius r. 

Also shown in the bottom panel of Fig. [2] is the local 
convective turnover time defined as 


_ 

'l~conv — ? 

t^conv 


(4) 


where Hp is the pressure scale height, and Vconv is convective 
velocity obtained from the mixing-length theory. Note that 
Tconv decreases toward the inner boundary at logT « 6.2; 
this is because the fractional radius of the convective inner 
boundary is very small (r/i? « 0.006, while Mr/M « 0.28). 
Fig. [2] indicates that in a considerable fraction of the zone 
where the radial displacement has a large amplitude, Tconv 
is much larger than the period and most of the energy is 
carried by radiation (i.e., Lrad/Atot ~ 1). This means that 
our approximation neglecting the divergence of convective 
flux is likely to be reasonably good. 

The absolute value of the imaginary part of the non- 
dimensional eigenfrequency w is comparable with the real 
part, indicating that the growth time is comparable with 
the pulsation period. In other words, the oscillatory con¬ 
vective modes are violently excited so that they should be 
capable of causing semi -regular light/velo c ity variations of 
the s ize observed (e.g., IWood et al.l Il999l : ISoszvhski et al.l 
l2004h during LSP variations. 

Fig.[3]shows the real part (a;,-; upper panel) and imagi- 
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nary part {uji ; lower panel) of normalized eigenfrequencies of 
oscillatory convective modes as a function of the luminosity 
of AGB models of various masses and with mixing-length 
parameters a = 1.2 (solid lines) and 1.9 (dashed lines). (We 
have selected only non-flashing, i.e., slowly evolving, AGB 
models.) Generally, ujr increases and \uji\ decreases with lu¬ 
minosity, because non-adiabaticity is stronger (or thermal 
time is shorter) in more luminous models. Conversely, as 
luminosity decreases, uir decreases and become essentially 
equal to zero at a certain luminosity for a given mass and 
a, below which the mode behaves as a monotonic g~ mode. 
This means that there is a minimum luminosity for the pres¬ 
ence of oscillatory convective modes depending on the mass 
and a. The minimum luminosity is higher for larger mass 
for a given a, while it is higher for larger a for a given 
mass (Fig. [3]). These dependences can be understood using 
the thermal time given in eq. m which should be inversely 
proportional to the strength of nonadiabaticity. For a larger 
mass model, a higher luminosity is needed to compensate a 
larger numerator, M — Mr, of eq. ( 0 ); the same is true for 
a model with larger a, because radius is smaller and hence 
mean density of the envelope is larger, which means that 
M — Mr is larger for a given temperature. 

As seen in Fig. [31 i^r increases as luminosity increases 
and become nearly equal to the absolute value of the imag¬ 
inary part, \u}i\, where each line is terminated. We have ter¬ 
minated the sequence because the amplitude in the radiative 
core becomes extremely large (> 10 ^) and rapidly oscillating 
for a mode with ujr > \^i\, although the amplitude distribu¬ 
tion in the envelope hardly differs from that shown in Fig. [31 
Such a huge amplitude in the core is caused by a coupling 
with a high-order g-mode; i.e., the mode is a mixed mode 
consisting of an oscillatory convective mode in the outer 
envelope and high-order g mode in the core. It is not pos¬ 
sible to resolve well such a large amplitude eigenfunction in 
the core (where r/R < 10~^ ) with a few ten thousands 
mass shells we are using. In such a mode, it is expected 
that non-linear effects in the core would become significant 
before the surface amplitude becomes appreciable. For this 
reason we consider only pure oscillatory convective modes 
with uJr < |<^i| to be observable. 

Nonradial pulsations in a red giant are known to 
have mixed mode properties, with p-mode characteristics 
in the envelope and g-m ode characteristics in the core 
dDziembowski et al.ll200ll ~l. The trapping of the p-mode in 
the envelope becomes very efficient as the luminosity in¬ 
creases as discussed bv lDupret et al.l (l2009l l. However, the 
coupling between an oscillatory convective mode in the en¬ 
velope and a core g-mode is different from the p- and g-mode 
coupling. In the case of the oscillatory convective mode, no 
isolation of the envelope oscillation from the core g-mode 
mode occurs if ujr > i.e., the amplitude of the coupled 

g-mode in the core is always exceedingly large. The differ¬ 
ence comes from the absence of an evanescent zone in the 
case of an oscillatory convection mode, while in the p- and 
g-mode coupling there is a narrow evanescent zone that sep¬ 
arates the p-mode propagating envelope from the g-mode 
propagating core. We note that the isolation of the oscilla¬ 
tory convective mode in the envelope does occur in helium 
shell flashing (thermal pulsing) models, in which a shell con¬ 
vection zone is present associated with the helium burning 
shell. The shell convection zone prevents the g-mode oscilla¬ 


tion from propaga ting into the core just as in the blue s uper- 
giants SPB stars (ISaio et al.ll2006l : Godart et al.ll2009l l. We 
have not included such flashing models in this paper be¬ 
cause the short time duration of shell flashing means such 
stars will be rare. 


4 COMPARISON WITH OBSERVED LSP 

VARIABILITY 

4.1 Period-luminosity relation 

Fig. [D compares period luminosity (PL) relations of red gi¬ 
ants in the LMC with theoretical period-luminosity relations 
of oscillatory convection modes (circles) and radial pulsa¬ 
tions (solid and dashed lines). Models with a = 1.9 are 
shown in the left panel and models with a = 1.2 are shown 
in the right panel. Our aim here is to match the periods of 
oscillatory convection modes with those of the LSPs belong¬ 
ing to sequence D while simultaneously matching the radial 
pulsation modes to sequence C (fundamental radial mode) 
and the shorter - perio d pulsation sequences. Period data 
(ISoszvhski et al.l l2009l 'l for red-giant variables were taken 
from the OGLE-III website (http://ogle.astrouw.edu.pl/). 
The luminosity of each star was obtained in the same way 
as for the sequence D stars shown in Fig. [T] 

The pulsation period is proportional to oc fiJr 
(eq. [T]). The PL relations for radial pulsations are roughly 
straight lines because uJr of each mode does not vary much 
along the AGB evolution. In contrast, the PL relation of the 
oscillatory convective mode has a peculiar shape. It bends 
and becomes nearly horizontal in less luminous parts, having 
a minimum period and a minimum luminosity. (The maxi¬ 
mum luminosity corresponds to the point where ujr = \u}i\ 
as discussed in the previous section.) The bending in the PL 
relation is caused by a rapid change in uJr as a function of lu¬ 
minosity near the minimum luminosity (see Fig. 0). Around 
the minimum luminosity, small Ur governs the PL relation; 
period decreases rapidly as luminosity increases due to an 
increase in At a certain luminosity, however, the R^'^ ef¬ 
fect exceeds the uJr effect so that the period starts to increase 
with luminosity (R). 

For models with the mixing-length parameter a = 1.2 
(right panel of Fig. [Tjl, the PL relation of the oscillatory 
convective modes above the ‘bend’ is consistent with se¬ 
quence D i.e., periods, the gradient of the relation, and the 
lower luminosity bound of sequence D agree with models of 
M > 1.0 Mq. However, these models are cooler than the 
red-giant sequence of the LMG in the HR diagram as seen 
in Fig. [T] and the radial mode periods, at least for masses 
less than ~ 2 M 0 , are longer than the observed periods. 

For models with a = 1.9 (left panel of Fig. [4|), which 
best agree with the distribution of red giants on the HR di¬ 
agram (Fig. [T]) in the LMC and which provide a reasonable 
match between radial-mode periods and observed pulsation 
periods, the periods of oscillatory convective modes above 
the ‘bend’ are too short for sequence D by a factor of ~ 2. Pe¬ 
riods can cross sequence D only below the ‘bend’, for which 
LUr < O.SoJi (open circles). Also, models of M ^ 1.0 Mq can¬ 
not explain the lowest part of the sequence D; the minimum 
luminosity for 1.0 Mq is logL/L 0 « 3.3, which is higher 
than the lower bound of the sequence D (log L/Lq « 2.8). 
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Figure 4. Period-luminosity relations of dipole oscillatory convective modes (circles) and radial modes in AGB models (non-flashing 
phases) with mixing length parameters of 1.9 (left) and 1.2 (right). Filled circles are for the modes with dUr > 0.5|a;i| or vise versa for open 
circles. Solid lines, long dashed lines, and short dashed lines correspond to the radial fundamental, first overtone, and second overtone 
modes, respectively. A maximum luminosity of each period-luminosity relation for dipole oscillatory convective modes arises because in 
more luminous models coupling with a high-order g-mode in the core is significant and hence the amplitude in the core is extremely 
high. Also shown are red-giant variables in the LMC with LSPs on sequence D (small black dots) and pulsating re d-giant variables with 
periods on sequences C, C', B and A (small green dots). The periods come from OGLE III llSoszvhski et aUl2009l l. 


Although the minimum luminosity for ~ 0.8 M© models de¬ 
creases as low as the lower bound of the sequence D, the 
ages of these models exceed the cosmic age. 

The problems in the a = 1.9 models are caused, at 
least partially, by the weak nonadiabaticity. The strength 
of nonadiabatic effects in a convection zone depends on the 
distribution of density and temperature, which in turn de¬ 
pends on the convection model. In this paper, we have used 
the mixing-length theory, which is a poor approximation 
for the convection in luminous red-giants, and we have ne¬ 
glected turbulent pressure and possible overshooting. A bet¬ 
ter convection theory including turbulent pressure and pos¬ 
sible overshooting might increase nonadiabatic effects sig¬ 
nificantly in outer layers, which could extend the luminosity 
range in which the oscillatory convective modes are found. 

4.2 Phase relations 

Fig. 0 shows predicted variations in luminosity, and pho- 
tospheric density (top panels), temperature and radial dis¬ 
placements (bottom panels) as functions of pulsation phase 
of dipole (t = 1) modes for three 2 M© models with differ¬ 
ent equilibrium luminosities; around the minimum period, 
middle, and the top of the PL sequence. The photospheric 
variations of 5T/T, 5p/p and SR/R were calculated using 
the eigenfunction at the mass shell corresponding to the 
photosphere in the equilibrium model, while 5L/L stands 
for the luminosity variation at the outer boundary which is 
located at an optical depth of 10“"^. We note that the quan¬ 
tities shown in Fig. [5] are not averaged across the surface 


but are local values proportional to the spherical harmonic 
17^(0, (?i). In particular, a dipole mode does not change the 
spherical shape and the mean radius of the stellar surface, 
although the local distance from the center of mass varies 
during pulsation. For this reason, we do not relate the ob¬ 
served radial velocity variation with a variation in the mean 
radius. 

Fig. [5] shows that temperature is higher when density 
perturbation is negative (positive buoyancy) and radial dis¬ 
placement is positive for oscillatory convection modes. These 
phase relations are similar to those for normal convective 
instability, but are quite different from ordinary p- and g- 
modes. Because of the peculiar phase relations in oscillatory 
convective modes, temperature and luminosity are higher 
when radial displacement is positive. The amplitude of 5TjT 
relative to the radial displacement decreases with luminos¬ 
ity; the ratio is about 1.5 in a model at logL/L© = 3.26 
(around the minimum period), while the ratio is reduced to 
about 0.2 at logL/L© = 3.78 (around the top of the PL 
relation). Non-adiabatic effects (thermal diffusion), which 
are stronger in more luminous models, reduce temperature 
variations. In addition, because of the stronger nonadiabatic 
effect, the phase relation between 5T/T and 5p/p shifts and 
the growth rate rj decreases as the luminosity increases. 

Fig. E] shows systematic trends of the amplitude and 
phase of STjT and SL/L as a function of luminosity for 
models of mass 2.0 (solid line), 1.3 (dashed line) and 1.0 M© 
(dot-dashed line). The behaviour relative to the luminosity 
at minimum period (filled circles) is similar among the cases 
with different masses. For a given mass, as the luminosity 






6 


H. Saio, P.R. Wood, M. Takayama, and Y. Ita 



2 r I I I I I I I I I I I I I I I I I I I 

0 0.2 0.4 0.6 0.8 1 


Phase 



0 0.2 0.4 0.6 0.0 1 

Phase 



2 r I I I I I I I I I I I I I I I I I I I 

0 0.2 0.4 0.6 0.8 1 

Phase 


Figure 5. One cycle of the photospheric temperature variations,(ST/T, radial displacements, SR/R (bottom panels), and photospheric 
density variation, Spjp, and luminosity variation at the outer boundary, SLjL (top panels) as a function of pulsation phase for dipole 
oscillatory convection modes in three models of 2.0 Mq {a = 1.2) with different luminosities, where n(d) is period in days, and r] is the 
growth rate defined as —LOiluir- The left panel is for a model around the minimum period (Fig.|4]|, the middle panel is for the same model 
as in Fig. [21 and the right panel is for a model close to the top of the period-luminosity relation in Fig. [J] The amplitude is normalizes 
as SR/R = 1 at zero-phase, and is proportional to a spherical harmonic Y-//^(9,(f>) with 9 and 0 being co-latitude and azimuthal angle. 


increases, the amplitudes of ST/T and SL/L decrease (rel¬ 
ative to the amplitude of SR/R). The phase of ST/T in¬ 
creases gradually with luminosity, while the phase of SL/L 
stays very small (i.e. SL is in phase with SR). The phase 
difference between ST/T (at photosphere) and SL/L (at the 
outer boundary) arises because of the entropy variation be¬ 
tween the photosphere and the outer boundary. The phase 
of the temperature variation at the outer boundary is very 
clos e to that of SL /R _ 

iNicholls et al.l (l2009l ) found various statistical proper¬ 
ties of the LSP variability. They found that radius and tem¬ 
perature variations based on the radial velocity variations 
are very different from those determined from photometry 
or spectroscopic analysis (which would be expected if the 
pulsations are dipole modes). Here we summarize the prop¬ 
erties of the majority of their sample based on photometric 
and spectroscopic analyses (i.e. we do not consider velocity 
data) as follows: (1) the phase delay of minimum light from 
the phase of minimum radius is mostly less than 0.2 period, 
(2) ATeff/Tefi < 0.03 and AR/R ~ 0.05. Q 

In the oscillatory convective modes shown in Fig. O 
minimum light occurs around minimum displacement which 
agrees reasonably well with the property (1). Because the 
pulsation is extremely nonadiabatic, temperature variation 
is very small. The ratio of the amplitudes of temperature and 
radius variations in the oscillatory convective mode ranges 
from about 2 to 0.2 ( from around minimum period to the 
top of the sequence; see Figs. [SJandO, while the property 
(2) above indicates the ratio is < 0.6 roughly consistent 
with the theoretical prediction. We note that if we use only 
photometric data, which cover all pulsation phases, the av¬ 
erage of the whole sample is found around 0.23. Thus, the 
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Figure 6. Amplitude (upper panel) and phase (lower panel) 
for the luminosity at the outer boundary (blue) and the pho¬ 
tospheric temperature (red) variations of oscillatory convective 
modes {£ = 1) versus luminosity. Solid, dashed, and dot-dashed 
lines are for 2.0,1.3, and 1.0 Mq models, respectively. The am¬ 
plitudes are normalized by the radial displacement, 5R/R. The 
phase is in units of period, with zero phase corresponding to the 
phase of maximum SR/R. A filled circle indicates the position of 
the minimum period for each case. 


^ We note that although the whole surface averages of ATeff and 
AR/R should be zero for nonradial pulsations, observational val¬ 
ues should be non-zero because they correspond to the means 
across the visible hemisphere; the cancellation effect is smallest 
in dipole modes among nonradial pulsations. 


properties of oscillatory convective modes seem consistent 
with the properties of the majority of the LSP sample of 

INicholls et al.l (l2009l). _ 

Recently, Takavama. Wood &: Ital (|2015l) examined 
color-magnitude variations during LSP variations for some 
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luminous SMC red giants. They found that matter should 
be ejected near the beginning of light decline. The mat¬ 
ter mi ght be related to the presence of chromosphere ac¬ 
tivity (IWood. Olivier fc Kawal^ |2004| ) and might be the 
origin of th e excess mid-IR emission observed around stars 
with LSPs l|Wood fc Nichollsl l2009ll . Fig. [5] indicates that 
maximum light occurs around maximum radial displace¬ 
ment, where mass ejection might be expected to occur. 
This is in agreement with phase of ejection found by 
iTakavama. Wood &: Ital ll2015l ~l. 


5 CONCLUSION 

We discussed properties of the oscillatory convective modes 
in luminous red giants and compared them with the prop¬ 
erties of the LSPs of variable red giants on the sequence D 
in the LMC. 

The phase relations among temperature, radial dis¬ 
placement, and density variations in the oscillatory convec¬ 
tive mode are similar to the relation expected for convective 
eddies rather than those of ordinary p- or g-modes; i.e., a 
positive temperature variation is associated with positive 
radial displacement and negative density variation. In addi¬ 
tion, because of strong nonadiabatic effects the amplitude 
of the surface temperature variations is reduced relative to 
the radial displacement. These properties of the oscillatory 
convective mode roughly agree with the properties of the 
LSP variations of red giants in the LMC. 

The period-luminosity relation of the oscillatory con¬ 
vective modes is found to be consistent with the observed 
relation of the LSPs of red-giants along sequence D if we 
adopt a mixing-length as small as 1.2 times pressure scale- 
height {a = 1.2). However, these models are much cooler 
than the AGB stars in the LMC. For models with a = 1.9, 
which agree with the effective temperatures of red giants 
in the LMC, periods are too short and the lower bound 
of the luminosity is too high, inconsistent with sequence 
D. Relatively weak nonadiabatic effect in the models with 
a = 1.9 would cause, at least partially, the discrepancy. The 
shortcomings of the present models might be improved by 
employing a more sophisticated convection theory including 
turbulent pressure and possible overshooting. We leave the 
problem for future investigations. 
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